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Abstract 

We present a search for the decays — > TT^i~^i~ and B^ — > TT^i~^£~, where £~^i~ is either a 
fj-~^fM~ or e^e~ pair, with a data sample of 657 mihion BB pairs collected with the Belle detector 
at the KEKB e^e~ collider. Signal events are reconstructed from a charged or a neutral pion 
candidate and a pair of oppositely charged electrons or muons. No significant signal is observed 
and we set an upper limit on the isospin-averaged branching fraction B{B — > 7r£~^£~) < 6.2 x 10~® 
at the 90% confidence level. 



PACS numbers: 13.25 Hw, 13.20 He 



In the Standard Model (SM), the decays 
B — s> TTi^i~ are forbidden at tree level and 
only proceed through flavor- changing neutral 
currents (FCNC). The possible lowest-order 
b — s> di^i~ processes are either a Z/7 pen- 
guin or a box diagram, as shown in Fig. [TJ 
FCNC b (s, d)i~^i~ decays can provide 
stringent tests of the SM in the flavor physics 
sector. For example, in the future it may 
be possible to measure the dilepton invari- 
ant mass spectrum and the forward-backward 
asymmetry in the center-of-mass system of 
the lepton pair, which can provide informa- 
tion on the coefficients of new operators as- 
sociated with theoretical models Most 



of the experimental studies and theoretical 
predictions are focused on b ^ s£+£~ de- 
cays jl. 111; however, signatures due to new 
physics may be observed in 6 — > di^i~ de- 
cays even if no such signature is found in 
b s£+t. 

The process b — > d£^i~ is suppressed 
by a factor of |Vjd/Vjsp ~ 0.04 relative to 
b s£~^£~, where Vtd and Vts are the el- 
ements of the Cabibbo-Kobayashi-Maskawa 
quark- mixing matrix [3]. The correspond- 
ing branching fraction of the exclusive — 
'7i~^i~^i~ mode is expected to be 3.3 x 10~^, 
which agrees with the simple power counting 
expectation . Assuming isospin symmetry. 
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the branching fraction for Ti^i^i~ is 

smaller by another factor of 1/2 x r^o / tb+ 
for these decays, where r^o and tb+ are the 
neutral and charged 5-meson lifetimes, re- 
spectively. Thus for SM branching fractions, 
only a few signal candidates are expected in 
our data sample. 

In this paper, we perform a search for the 
B+ n+i+i- and 5° ^ n^i+fr decays, 
where stands for a ^~^fi~ or e"'"e~ pair. 
A data sample of 657 million BB pairs col- 
lected with the Belle detector at the KEKB 
e~^e~ collider is examined. Charge-conjugate 
decays are implied throughout this paper. 




FIG. 1: The quark-level diagrams for B — > 
7ri~^£~ decays. 

The Belle detector is a large-solid-angle 
magnetic spectrometer located at the KEKB 
collider , and consists of a silicon vertex de- 
tector (SVD), a 50-layer central drift cham- 
ber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like ar- 
rangement of time-of-fiight scintillation coun- 
ters (TOF), and an electromagnetic calorime- 
ter comprised of CsI(Tl) crystals (ECL) lo- 
cated inside a superconducting solenoid that 
provides a 1.5 T magnetic field. An iron 
flux-return located outside the coil is instru- 
mented to detect mesons and to identify 
muons (KLM). The detector is described in 
detail elsewhere Ml. 



We consider all tracks that have a maxi- 
mum distance to the interaction point (IP) of 
5 cm in the beam direction (z) and of 0.5 cm 
in the transverse plane {r-cj)). We select vr^ 
candidates from charged tracks having a kaon 
likelihood ratio less than 0.4; the kaon likeli- 
hood ratio is deflned by TZk = C-k/ {Ck+C.„), 
where Ck (jC^t) denotes a likelihood that com- 
bines measurements from the ACC, the TOF, 
and dE/dx from the CDC for the (tt^) 
tracks. The selection efficiency is about 89% 
for pions while it removes 91% of kaons. In 
addition to the information included in the 
kaon likelihood ratio, muon (electron) candi- 
dates are required to be associated with KLM 
detector hits (ECL calorimeter showers). We 
define the likelihood ratio TZx {x denotes ix 
or e) as = C^/{C^ + £„ot-x), where £^ 
and Cnot-x are the hkelihood measurements 
from the relevant detectors [8|. We select 
ji^ candidates with IZ^ > 0.9 if the momen- 
tum of /i^ is greater than 1 GeV/c; for the 
/i^ candidates with lower momentum (0.7- 
1.0 GeV/c), IZ^ is required to be greater than 
0.97. These requirements retain about 80% of 
muons while removing 98.5% of pions. Can- 
didates for are required to have a mini- 
mum momentum of 0.4 GeV/c, TZ^ > 0.9, and 
TZ^ < 0.8. These requirements keep about 
90% of electrons while they remove 99.7% of 
pions. Bremsstrahlung photons emitted by 
the electrons are recovered by adding neutral 
clusters found within a 50 mrad cone along 
the electron direction. The energy of the ad- 
ditional photon is required to be less than 
0.5 GeV. For 7r° 77 candidates, a mini- 
mum photon energy of 50 MeV is required 
and the invariant mass must be in the range 
115 MeV/c^ < M(77) < 152 MeV/c^, corre- 
sponding to ± 3 (T for the reconstruction. 
Requirements on the photon energy asymme- 
try, \E] - E^\/{E] + E^) < 0.9, and the min- 
imum momentum of the 7r° candidate, p^o > 
200 MeV/c, are introduced to suppress the 
combinatorial background. 

5-meson candidates are reconstructed in 
one of the following modes: B~^ 7i~^e'^e~, 
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B+ 7r+/i+/i-, 5° ^ 7r°e+e-, or 5° ^ 
7r°/i+/i~. Signal candidates are selected 
using the beam-energy constrained mass 
^bc = \/ -^beam - p| ^nd the energy differ- 
ence l\E = Eb — -Eboam, where Eb and 
are the reconstructed energy and momentum 
of the B candidate in the T(45') center-of- 
mass (CM) frame, and -Ebeam is the beam 
energy in this frame. We require S- meson 
candidates to be within the region Mbc > 
5.20 GeV/c2 and -0.1 GeV < < 0.3 GeV 
(-0.15 GeV < < 0.3 GeV) for the B+ ^ 
7r+£+£- (5° ^ 7r°£+£-) decays. The signal 
region is defined by 5.27 GeV/c^ < A/bc < 
5.29 GeV/c2, -0.035 (-0.08) GeV < AE < 
0.035 GeV for the 7r+yU+/i~ (n^fi'^ii^) mode, 
and -0.055 (-0.1) GeV < AE < 0.035 
GeV for the vr+e'^e" {ir^e^e") mode. The 
lower AE bound for the candidate region is 
designed to exclude possible contamination 
from similar B decays, e.g. B K*ii, while 
the wider region in AE for B^ n^i^i' de- 
cays is chosen to include the long tail of the 
signal. 

The dominant source of background is 
continuum e^e~ — > qq events (g = u, d, c, s). 
A Fisher discriminant including 16 modified 
Fox- Wolfram moments is used to exploit 
the differences between the event shapes for 
continuum qq production (jet-like) and for 
B decay (spherical) in the e~^e~ rest frame. 
We combine 1) the Fisher discriminant, 2) 
the angle between the momentum vector of 
reconstructed B candidate and beam direc- 
tion {cos6b), and 3) the distance in the z 
direction between the candidate B vertex 
and a vertex position formed by the charged 
tracks that are not associated with the can- 
didate 5-meson into a single likelihood ra- 
tio TZ = Cs/i^s + ^qq), where Cs i'^qq) de- 
notes the signal (continuum) likelihood. The 
background from B — > iuX decays is also 
found to be large. An additional likelihood 
ratio TZb = £s/(£s + C^-^) including cos 9b 
and the overall missing energy is introduced, 
where C^-j^ is the likelihood for BB events. 

Continuum background suppression is im- 



proved by including S-flavor tagging infor- 
mation [lO|, which is parameterized by a dis- 
crete variable gtag indicating the flavor of the 
tagging _B- meson candidate and a quality pa- 
rameter r (ranging from for no flavor in- 
formation to 1 for unambiguous flavor as- 
signment). Selection criteria for TZ and TZb 
are determined by maximizing the value of 
S/ VS + B, where S and B denote the ex- 
pected yields of signal and background events 
in the signal region, respectively, in different 
^rec ■ qtag ■ f regions, where q^ec is the charge 
of reconstructed B candidate. Events with 
qrcc ■ qtag ■ r close to —1 are considered to be 
well-tagged and are unlikely to be from con- 
tinuum processes. For neutral 5-meson de- 
cays, only the dependence on r is considered. 

The decays B J/ipX and ip'X are the 
dominant peaking backgrounds in the AE- 
Mbc candidate region. Events in the regions 

- 0.10 < [M(/i+/i-) - m{J/tp, tp')] < 0.08 , 
-0.20 < [M(e+e-) - m{J/ip, < 0.07 , 

0.86 ■ [M (/i+/i-) - m{J/ip, ip') - 0.08] < 

AE/c^ < 0.86 ■ [M(/i+/i-) - m{J/ij, i^') + 0.10] , 

0.94 ■ [M{e+e') - m{J/^, - 0.07] < 

AE/c^ < 0.94 ■ [M(e+e-) - m{J/ip,ip') + 0.20] . 

in the AE-Mi+i- plane are rejected (the 
limits on the regions are given in units of 
GeV/c2, see Fig. Wj- The decay 5+ ^ 
J/ijj{ip')h^ {h^ = K^,7T~^) can also con- 
tribute to the B^ 'K~^fi~^fi~ sample. If 
a muon from J/tpltp') is misidentified as a 
pion and another non-muon track is at the 
same time misidentified as a muon, such 
a background event cannot be removed by 
the criteria described above. Thus, in addi- 
tion, events with -0.10 GeV/c^ < M(/i7r) - 
m{J / Ip ^ Ip') < 0.08 GeV/c^ are removed from 
the B^ — ^ ir'^fi'^fi" sample. An event re- 
constructed B^ D'^{^ TT+TT )7r+ or 
B^ — * /}"(—> 7r°7r~)7r+ decay can contribute 
to the B —>■ 7r/i^yU~ sample if both pions 
are misidentified as muons. We reject these 
events by requiring |M(7r/i) — m{D)\ > 0.02 
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GeV/c^. The invariant mass of the electron 
pair is also required to be greater than 0.14 
GeV/c^ in order to remove background from 
photon conversions and from vr^ 7e^e~ 
decays. The background contributions from 
B J/'ip{'ip')X and other b ^ c decays 
are estimated using large Monte Carlo (MC) 
samples. Background from charmless three- 
body B decays, such as B^ Tr+Tr^Tr"*" and 
B^ vr^Tr'Tr"*" is estimated using the mea- 
sured data, taking into account the probabili- 
ties of the pions being misidentified as muons. 
The results are consistent with MC simula- 
tions. The background yield from misidenti- 
fied B~^ — >■ K~^i~^i~ decays is calculated from 
the yield of a dedicated B~^ K~^i~^i~ anal- 
ysis with a similar analysis procedure, which 
is then scaled by the kaon to pion misidenti- 
fication rate. 




M^+^- (GeV/cl 




2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 

M3+ - (GeV/c') 



FIG. 2: AE vs. M^+^- distribution of B ^ 
J/tpX and ^l^'X MC events feeding into Tr^fi~^fi~ 
(left) and Tr^e'^e~ (right) modes, respectively. 

If there are multiple candidates in an 



event, the candidate with the smallest vertex 
fit is chosen for the B^ — s> 7r~^i~^i~ modes, 
and the one with the best TZ value is selected 
for the B^ iT^i'^i~ modes. The fraction 
of such events are about 11%-13% and 16%- 
20% for B+ 7i+i+i- and B^ n^f+i-, 
respectively, according to a MC study. 

We perform a simultaneous extended un- 
binned maximum likelihood fit to AE and 
Mbc with the following likelihood function: 



C 



-(Ns+Nqq+N^-X+NKU + Nhnn) 



N\ 



N 



NbbPbb + 



i=l 



NccxPccX + ^KuPru + A^'hTTTT-PLTr] • 

where denotes the number of observed 
events in the candidate region, and A^^ 
{PI), N,-, (P^,), iV^s (P^^), iV,,x {P:-.x). 
Nku {Pru)^ and Nh^^ (PLn) denote the 
event yields (the probability density func- 
tions, PDFs, for the i-th event) for signal, 
continuum, b ^ c decays, B J /ip{il}')X, 
background contributions from B K£~^i~ 
and B hn^TT' decays. Since contributions 
from B J/i){ip')X and B Ki+i- pro- 
cesses are found to be negligible for P° — > 
7r'^i~^£~ decays, the associated PDFs are ex- 
cluded from the fits. The signal PDFs for 
B~^ 7r^i~^i~ are assumed to be Gaussian 
in both AE and Mbc. The means and widths 
are verified using B —>■ J/ipK decays. Addi- 
tional Gaussian functions are used to model 
the tails in the AE distributions. For the 
P° TT^i+i- modes, the signal PDF is 
modeled by a two-dimensional smoothed his- 
togram function, including the dependence 
on M {£+£-). The continuum PDFs for Mbc 
and AE are represented by an empirical 
background function introduced by ARGUS 
[ll| and a second-order polynomial, respec- 
tively. The PDF for B J/i){ip')X decays 
consists of a peaking part and a non-peaking 
part. The peaking PDF is modeled by Gaus- 
sian functions in AE and Mbc. The non- 
peaking part, as well as the PDF for other 
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b ^ c decays, is described by the same PDF 
formulae used for continuum events, but with 
different parameter values. The distributions 
of AE and Mbc for the B Ki+i- and B 
hir^TT" contributions to the B'^ 7r~^i^i~ 
modes are found to be peaking and modeled 
by Gaussian functions. Due to a much larger 



tail on AE for the B" 



TT 



B ^ h 



TT ' TT 



contributions to B^ 



decays, the 



TT 



are modeled by smoothed histogram func- 
tions. 

Yields for signal and continuum, and the 
continuum PDF parameters are allowed to 
float in the fit while the yields and parameters 
for other components are fixed. The observed 
signal yields and branching fractions [B) are 
summarized in Table [H The significance is 
defined as -^Z— 21n {Co / £max) , where Cq is the 
likelihood with signal yield constrained to be 
zero, and £max is the maximum likelihood. 
The distributions of AE and Mbc with fit re- 
sults superimposed are shown in Fig. [3l 

Systematic uncertainties are summarized 
in Table HTl The dominant systematic uncer- 
tainties for the B~^ 7T^i~^i~ modes stem 
from electron (4.4%) and muon (5.5%) iden- 
tification efficiencies, uncertainties of MC de- 
cay models (2.6%-3.8%), and tracking effi- 
ciencies (3%). The signal MC samples are 
generated based on the decay model derived 
from [12], and the uncertainties are evalu- 
ated by comparing the MC and data yields 
for B+ K+i+i- function of dilepton 
invariant mass. Other uncertainties such as 
background suppression, which is studied us- 
ing B — ^ J/ipK decays, and PDF modeling 
are all small. For the B^ — > modes, 
the dominant uncertainties are from lepton 
identifications, PDF modeling (5.3%), vr^ re- 
construction efficiency (4%), and MC models. 
Other systematic uncertainties are found to 
be small. There are also some small additive 
uncertainties that are not proportional to the 
signal branching fractions. An example is the 
uncertainty in the background level, which is 
estimated to be about 0.3-0.8 events. In ad- 
dition, a fitting bias is found in the MC sim- 



ulation when signal yields are small (fewer 
than 2 events). Thus, we quote an additional 
uncertainty of 0.5 events for each mode. 

As no statistically significant excess of sig- 
nal events is observed, we calculate upper 
limits of the branching fractions by integrat- 
ing the likelihood function up to the 90% con- 
fidence level. We find. 



B{B+ - 
B{B+ - 
BiB"^ - 
BiB"" 



< 6.9 X 10"^ , 

< 8.0 X 10"® , 

< 18.4 X 10"® , 

< 22.7 X 10"® . (1) 



Systematic uncertainties are included by con- 
voluting the likelihood function with a Gaus- 
sian function. The lepton-fiavor combined 
limits are given by 

B{B+ Tx+tr) < 4.9 X 10"® , 
i3(5° ^ TxH+t) < 15.4 X 10"® . (2) 

Assuming isospin symmetry, the average 
limit for B 7r£+£~ is 



B{B^ 



TT 



') < 6.2 X 10^ 



In conclusion, the upper limit on the 
isospin-averaged branching fraction is about 
twice the SM expectation, and two thirds of 
the BaBar measurement 13|. A much larger 
data set, such as will be available from the 
proposed super B factory [3], is needed to 
observe this decay mode if the branching frac- 
tion is at the level predicted by the SM. 

We thank the KEKB group for the excel- 
lent operation of the accelerator, the KEK 
cryogenics group for the efficient operation 
of the solenoid, and the KEK computer 
group and the National Institute of Infor- 
matics for valuable computing and SINET3 
network support. We acknowledge support 
from the Ministry of Education, Culture, 
Sports, Science, and Technology of Japan 
and the Japan Society for the Promotion 
of Science; the Australian Research Council 
and the Australian Department of Education, 
Science and Training; the National Natural 
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FIG. 3: Distributions of AE (Mbc) with fit results superimposed for the events in the Mbc (AS) 
signal region. The solid curves represent the fit results, while the solid (dashed) peaks represent 
the signal {K^i~^i~ background) component. 



TABLE I: A summary of the signal yields (Ng), reconstruction efficiencies (e), statistical significance 
(S), branching fractions (B), and the corresponding upper limits (U.L.) at the 90% confidence level. 
The errors associated with branching fractions are sequentially statistical errors, systematic errors 
proportional to the branching fractions (given in Table II), and additive systematic errors related 
to yield extraction. 



Mode 




Ns 


e(%) 


S 


B (10-8) 


U.L. (10-S) 


B+- 






13.1 


0.2a 


0.6l|^±0.0±0.7 


6.9 


B+ - 






13.8 


0.6a 


1.5l:|^±0.1±0.8 


8.0 


B+ - 








0.6a 


i.ol:^J±o.i±o.4 


4.9 


i?o^ 




r: 1+4.2 
'-'•-'--3.3 


9.6 


1.8a 


8.ll|^±0.8±1.0 


18.4 


B^^ 




9 7+5.2 
^- '-4.0 


7.4 


0.7a 


b.f)tf£±0.b±2.0 


22.7 


B^^ 








2.0a 


7.4l|^±0.7±0.8 


15.4 


B Tri+£- 






1.2a 


2.4+|^±0.2±0.2 


6.2 
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TABLE II: Summary of the contributions of systematic uncertainty proportional to the branching 
fractions (in %). 



Decay 


TT+fX+II- B + 


TT+e+e" B^ - 


Tr° l_r B^ - 


n _i_ 

7r°e+e" 


Tracking efficiency 


3.0 


3.0 


2.0 


2.0 


e / //-identification 


5.5 


4.4 


5.5 


4.4 


7r'''-id./7r'' reconstruction 


1.0 


1.0 


4.0 


4.0 


Background suppression 


1.8 


1.8 


1.8 


2.4 


PDF Modefing 


2.2 


2.2 


5.3 


5.3 


MC decay model 


3.8 


2.6 


3.8 


2.6 


MC statistics 


0.8 


0.8 


1.1 


1.1 


N{BB) pairs 


1.3 


1.3 


1.3 


1.3 


Total 


8.1 


6.8 


10.0 


9.1 



Education of Taiwan; and the U.S. Depart- 
ment of Energy. 
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